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Introduction 

The flow over the two-dimensional hump model is computed by solving the RANS equations 
with k-co (SST) model. 

Solution Methodology 

The governing equations, the flow equations and the turbulent equations, are solved using the 5th 
order accurate weighted essentially non-oscillatory (WENO) scheme for space discretization and 
using explicit third order total-variation-diminishing (TVD) Runge-Kutta scheme for time 
integration. The WENO and the TVD methods and the formulas are explained in [1] and the 
application of ENO method to N-S equations is given in [2], The solution method implemented 
in this computation is described in detail in [3]. 

Model used 

Standard k-co (SST) model is used and the equations and the model coefficients are described in 
[4,5,6], 

Implementation and Details 

The computational domain extends from x/c=-10. to 4. in the streamwise direction and extends 
from the splitter plate to the upper tunnel wall in the normal direction. The leading edge of the 
splitter plate is modeled as a super ellipse with 0.25 in. half thickness and an aspect ratio of 2. 
The leading edge is located at x/c=-5.9. C-Type grid is used around the splitter plate and a 
rectangular grid is added upstream of the leading edge as shown in Fig. 1. This grid overlaps the 
C grid and 5 th order central interpolation is used to transfer the flow variables from one grid to 
the other at the boundaries. (651*151) grid size is used around the splitter plate and the hump 
and (101*51) grid size is used in the rectangular region. 
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Figure 1: Overlapping grid near the leading edge of the splitter plate. 


Following boundary conditions are implemented at different boundaries: 
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At the upper wall inviscid conditions are applied. 
dp du dE 

dy dy dy 

At the lower wall viscous conditions are used. 
u = v = 0, 

T — T = T 

w free stream’ 
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and p is computed from the continuity equations. 

From the leading edge of the splitter plate to the inflow boundary symmetric conditions 
are used. 

At the inflow boundary stagnation pressure, one Riemann variable and normal velocity 
v=0 are prescribed and the second Riemann variable is solved for to obtain the other flow 
quantity. 

At the outflow boundary the pressure is specified to obtain the required Mach number 
and characteristic-type boundary conditions are implemented similar to as described in 
[7] to obtain other flow variables. 

In the suction case, boundary conditions are applied on the surface of the hump across the 
suction slot. The suction slot extends from x/c=.6541 to .6584 and across the slot normal 

mass flow rate is specified to the experimental value. A suction distribution of the form 
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is used. Other forms have been tried and all of them yield the same results for a fixed 
total suction rate. The other flow quantities are obtained using the characteristic type 
boundary conditions [7]. 


Following boundary conditions are implemented for the turbulent quantities at different 
boundaries. 

1. At the inflow boundary small values are prescribed for k and p T . 
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2. At the outflow boundary k and co are solved for from the governing equations. Higher 
order extrapolation condition is also tried and it gives the same results. 

3. At the lower viscous wall k=0 condition is used and following exact boundary condition 
is derived for co. 

4. In the suction case, across the suction slot linear extrapolation is used to obtain the 
turbulent quantities on the surface. 


Since the variable co becomes singular near a viscous wall, in practice a large approximate value 
is prescribed at the wall. 
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where d w is the distance to the first grid point from the wall. This is an approximate boundary 
condition and when it is implemented in the higher order scheme, oscillations and convergence 

problems are encountered and the following exact condition is derived for co wall . By realizing 
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type singularity for the variable co arises because of the balance between the 


dissipation term and the viscous diffusion term in the co equation, the singularity is removed by 
rewriting the variable co as 
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where y n is the normal distance to the wall, C is a constant and C0i is the new variable which is 
now regular near the wall. When this is substituted into the co equation the following equation is 
obtained for C0i, which is similar to the co equation except for the source term. 
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The value of C0i at the wall becomes 
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Here the variable co is nondimensionalised by 



( 8 ) 


(9) 


Hence the procedure is to use the C0i equation for the first few points near the wall and switch to 
the co equation away from the wall. In these computations (Oi equation is solved for the first ten 
points near the wall. Figure 2 shows the distribution of k, co and (Oi near the wall and it is seen 
that this technique resolves the viscous layer smoothly even though co is infinite at the wall. 
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Figure 2: Variation of k, co and C0i near the wall at x/c=0.4. 

Figure 3 shows the contours of the U velocity near the leading edge region and for the entire 
computational domain. Near the leading edge region the flow separates and forms a small 
separation bubble. In the no flow case, the solution converged for the flow equations and the 
turbulent equations very well. In the suction case, the maximum residual in the co equation near 
the suction slot edges converged only by three orders. This may be due to the boundary 
conditions used for the turbulent quantities near the suction slot. 
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Figure 3: Contours of U velocity near the leading edge and over the hump. 

References 

1. Shu, Chi-Wang, “ Essentially Non-Oscillatory and Weighted Essentially Non- Oscillatory 
Schemes for Hyperbolic Conservation Laws,” NASA/CR-97-206253 and ICASE Report 
NO. 97-65, 1997. 

2. Atkins, H. L., “High-Order ENO Methods for the Unsteady Compressible Navier-Stokes 
Equations,” AIAA Paper 91-1557, 1991. 

3. Balakumar, P., “Stability of Hypersonic Boundary-Layers Over a Compression Comer,” 
AIAA Paper 2002-2848, 2002. 

4. Menter, F. R., “Improved two-Equation k-co Turbulence Models for Aerodynamic 
Flows,” NASA Technical Memorandum, NO. 103975, 1992. 

5. Menter, F. R., “Zonal Two Equation k-co Turbulence Models for Aerodynamic Flows,” 
AIAA Paper 93-2906, 1993. 

6. Menter, F. R. and Rumsey, C. L., “Assessment of Two-equation Turbulence Models for 
transonic Flows,” AIAA Paper 94-2343, 1994. 

7. Sesterhenn, J., “A characteristic-type formulation of the Navier-Stokes equations for high 
order upwind schemes,” Computers & Fluids, 30, pp. 37-67, 2001. 


3.15.5 


